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Abstract 

The contribution of hadronic vacuum polarization to the hyperfine 
spUtting of the muonic hydrogen ground state is calculated with the 
account of experimental data on the cross section of e~^e~ annihila- 
tion into hadrons and the dipole parameterization of electromagnetic 
proton form factors. 
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Over many years the comparison of experimental and theoretical results 
for hyperfine splitting of hydrogen, muonium and positronium ground state 
energy level is considered as an important verification of our understanding 
of the bound state problem in quantum electrodynamics. Especially high 
precision was achieved at the present time when measuring hydrogen atom 
hyperfine structure [|l], 



Theoretical value of hyperfine splitting, which ordinary written in the form 



(/ip is the proton magnetic moment, expressed in nuclear magnetons; mp, 
rrie are proton and electron masses correspondingly), contains besides elec- 
tromagnetic correction 5^^^ some more two terms 5^ and 5^ . These cor- 
rections take into account the influence of strong interaction. The term 5^ 
describes the effects of proton finite size and recoil contributions, and 5^ is 
the correction of proton polarizability. Particle structure members 5^ and 
5^ are absent for purely leptonic systems 0. Electrodynamical correction 
^QED hyperfine structure of hydrogen atom, accounting relativistic and 
radiative effects, is known at present time with the accuracy 0,12 ppm [^]. 
The increase of accuracy in theoretical expression (2) is impossible without 
the detailed consideration of strong interaction in the terms 5^ and 5^ . The 
calculation of different contributions to the hydrogen atom hyperfine split- 
ting (HPS), connected with proton structure (correction 5'*)) was done in the 
paper [Q. Authors of have proposed an experimental method to obtain 
more exact information about proton electromagnetic structure. It based on 
the measurement of hyperfine splitting of the muonic hydrogen (yup) ground 
state with the accuracy 10~^. It is the authors opinion that such measure- 
ment will allow to find the relation between the corrections 5^ and 5^, what 
will increase the accuracy of theoretical formula (2). In the case of muonic 
hydrogen, when the light electron must be replaced by the heavy muon, the 
relative contribution of different effects in (2) is changed. So, the hadronic 
vacuum polarization contribution, which was omitted in (2) in view of it's 
small value by comparison with 5^ and 5^, may have more relative contribu- 
tion for muonic hydrogen HPS. The effect of hadronic vacuum polarization 
must be considered in order to extract correctly the parameters 5^ and 5^ 
from experiment. In this paper we have calculated the hadronic vacuum po- 
larization contribution to the hyperfine structure of muonic hydrogen on the 



AE^Yp' = 1420405.7517667(9)Mz. 



(1) 




(2) 
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basis of quasipotential method 0. In this approach the particle interaction 
operator for I7- and 27- processes takes the form [^: 

V = Vi^ + = V' + A\/, (3) 

Vi^ = Ti^, V2-y = — Ti^ X X Ti^, 

where V^"^ is Coulomb potential, Ti^, are one- and two-photon scattering 
amplitudes beyond energy surface, [G-^]"^ = —fP')/2nR is the free two- 
particle Green function {^r is relativistic reduced mass, 6^ is the square of 
particle relative motion momentum on energy surface). Then we can write 
the displacement of energy level as follows: 

=< ^^\/\V\^^ >, (4) 

where ip'^ is the ordinary Coulomb wave function. The hadronic vacuum po- 
larization contribution to the energy spectrum is determined by Feynman 
diagrams represented on Fig.l by means (3). To take into account vacuum 
polarization in these one-loop diagrams, we must do the following substitu- 
tion in photon propagator 

+ ie xn J Jsth. k"^ — s + ie^ 

where the spectral function p{s) is connected with known cross section of 
e~^e~ - annihilation to hadrons a^: 

, , R{s) a^ie^e" — ^ hadrons) , ^ 

3s 3scr^^(e+e ) 

and (Jy^^ie^e" — * /U^/x^) = 47rQ;^/3s is the e+e~- annihilation cross section to 
muonic pair. 

Let consider the proton factor, corresponding to two-photon exchange 
diagrams: 

^ r/-(7°^2-P + mp)r- r'-(7°E2+p + m,)r^ 

where j)^ is the loop 4-momentum, 4-momentum t = (0, p— q) is determined 
by relative motion momenta p, q of initial and final states correspondingly. 
Proton vertex function may be expressed through Dirac p(p^) and Pauli 
/(p2) form factors, describing the proton electromagnetic structure: 



(8) 
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Figure 1: Feynman diagrams, defining the hadronic vacuum polarization 
contribution to hyperfine structure of atom [fip) 

wliere p(0) = e, /(O) = en, k is tlie proton anomalous magnetic moment: 
fi;= 1.792847386 (63) |T|. Muonic tensor of considered one-loop amplitudes is 
equal: 

_ 7n7°^i+P + m,)7'- 
~ p2 + b^ + 2E^Po + ze' ^ ^ 

The hadronic vacuum polarization contribution to the ground state hyperfine 
splitting of muonium and positronium was obtained in 0]. It is clear that for 
exact calculation of discussed one-loop corrections in (/ip), we must consider 
the dependence of form factors p(p^), /(p^) from loop momentum p. To do 
numerical calculations of the corrections 0{a^) caused by hadronic vacuum 
polarization we can use dipole parameterization of proton electromagnetic 
form factors of the form: 

Ge - Gm - Ge 

PiP ) = — ^_ ^ — , fip ) = > (10) 

= Ge = —2, (11) 

_ £1] (1 _ £1 
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where the parameter of proton structure was taken in the kind A = 0, 898mp 
[||]. Multiplying the expressions (7) and (9) for construction of quasipotential, 
we obtain the following sum of form factor terms with the accuracy 0{a^): 



1 1 
+ 



rrir. 



D,{-p) D,{p)J D,{p) 



+4:p{p'^)f{p^) < pY' >< in^ > 



Dp{-p) Dp{p)J D^ip) 



+ 



rrir. 



p^ < ^^pY X YPl'^ > 



D,{-p) Dp{p)J D^ip) 



[12) 



2mpp < -f^T >< TPl > 



D,i-p) Dpip)J D^ip) 



+ (^^] 2my < YPl" >< YPl" > 
where 



rrir, 



Dpi-p) D,{p) D^{py 



Dp{±p) =p^ ± 2mpPo + ie, D^{p) = p^ + 2m^pQ + is. 

To make numerical calculations let transform (12) to convenient form. First 
of all we pick out in proton factor l/Dp{—p) of crossed two-photon diagram 
the separate term with 6{pq) function, which cancelled in future calculations 
by similar addendum of iteration part of quasipotential (3): 



Dp{-p) 2{mp + m 



-27ci6{po) 



D,{p) 



Po + ie {-po + ie)Dp{-p) 



■ (13) 



It is useful also to represent the other proton factor 1/Dp{p), appearing in 
(12), in the following manner: 



Dp{p) 2{mp - rUf,) {po + ie) 



1 



D,ip) 
DM 



(14) 



Accounting these relations, we may transform different parts of (12) as fol- 
lows: 



D,{p) 



Dp{-p) Dp{p) 



rrir 



ml -ml 



1 



[Dp{p) D,{p)\ 



(15) 
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ml -ml 



ml — ml 



2 

mt 



Dpip) 



m. 



(16) 



(17) 



After that we rotate the po contour to the imaginary axis and call the new 
variable ^. The integrated factors are modified as follows: 



1 



1 



1 



-p + ^ - 2imp^ 
fTie ~^+~^' Dp{p) ^ (p2 + ^2)2 + 4^V2- 



(18) 



As a result of such transformations we obtain that the contribution of 1, 3 
and 4 form factor terms in (12) to muonic hydrogen HFS with the account of 
hadronic vacuum polarization is determined by triple integral of the following 
form: 

32a2 



hf.s 



mpm^ 



mt 



ml vr^ 



SO 



p{s)ds / p dp 



(19) 



[p' - &y + erne + m + Afpf 



mt 



m. 



{p^ + ey + 4mie {p^ + + 4:ml^\ 



To do integration over s, we have used the parameterization of p{s) by means 
of new experimental data on annihilation cross section e~^e~ — > hadrons 
0' 11' 12' '^^^ main contribution to the cross section is determined 
by the process e+ + e~ 7r+ + tt^. The cross section of this reaction is 
proportional to the squared module of vr meson form factor F^^. Supposing 
that the p exchange diagram gives the basic contribution to F^, we have 
used the following well-known expression of pionic form factor, taking into 
account the contribution oi p — uj interference term fll 
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m^(l + dVp/mp 



[ml - s + impTp{p/po)^mp/ ^/s 



i(f> 



mr. 



mt 



s + im^T^ ' 



(20) 



where the p meson decay width rp=0.15 Gev, p meson mass mp=0.77 Gev 
and d=0.5. The effects of p — u) interference are described by the parameters 
from |TD[ . Substituting (20) to the spectral function 



(s) = 1- IF 



(21) 
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we obtain from (19) the corresponding contribution to HFS of muonic hy- 
drogen: 

A£;^-^i'(p ^ 27r) = 402.40 MHz, = 15796.2 GHz. (22) 

The term of (19), proportional to P /rn^ leads to the negative contribution 
23,84 MHz in the energy spectrum. It is interesting to compare the quantity 
(19) with the result of it calculation in the point-like proton approximation, 
when we assume p{k'^) ~ p(0) = e, /(fc^) ~ /(O) = en. Keeping in mind 
that the main contribution in (19) is determined by first and third terms 
of (12), let transform their quantity in the hyperfine splitting of (/xp), ful- 
filling analytical integration over loop momentum p by means of Feynman 
parameterization. As a result we find: 



-K J nip 



■Ef 



m. 



Ami 



dsp{s) / dx[S{mp,x) — S{mn,x)] , 



(23) 



where 



5'(m, x) 



16 — 6x — + 4k(8 — 6x + x^) 



, E^ = Sa^rnlrnp/{rnp + m,^y 



o I s(l — x) 



Integration over parameter x may be done analytically also, and to inte- 
grate over s, wc have used parameterization of p{s) in the form (21). Then 
the contribution of first and third terms of (12) to HFS in the point-like 
proton approximation is equal to: 



AEl^^-'(p 27r) = 882.95 MHz. 



(24) 



Comparing this value with the result (22), we may state that using of form 
factors (10), (11) in vertex operator (8) leads to essential expected decrease 
of corresponding correction by comparison to approximation p(p^) ~ p(0), 
/(p^) ^ /(O). The other parts of (12) give the contribution to HFS, which 
may be represented as a sum of two integrals similar to (19): 



9 99 
mi — mf. TT^ 



p{s)ds 



so 



I p'^dp 

1 



(25) 
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+ + s) 



1 



(p2 + ^2)2 + 4^2^2 (p2 + ^2)2 + 4^2^2^ 

We have represented in the table the total contribution of all form factor 
terms of (12) with function Pt^t^{^s) to the HFS of muonic hydrogen, which 
is formed by the expressions (19), (25) and (26). Let take into account the 
contributions of uj and mesons to hyperfine structure of (/^p). To do this 
would require Breit-Wigner representation of spectral function p(s): 

where = 0.782 Gev, = 1.019 Gev, = 8.43 Mev, = 4.43 Mev, 
= 0.60 Kev, = 1.37 Kev [|T^. Then after calculation of triple integrals 



(19) we obtain the summary contribution of the following kind: 

A£;^-''-"(cj) = 45.74 MHz, (28) 

AE^-^'(0) = 48.97 Mi/z. (29) 

Spectral density R(s) was parameterized in the energy range 0.64 < s < 
1.44, between uj and resonances on the basis of present experimental data 
[0 on e+e~ annihilation cross section in the form: R{s) = A + Bs + Cs^ 
(A=22.68, B=-37.72, 0=15.84). To obtain the contributions of J/^ and T 
particle families we have also used the Breit-Wigner formula (27) with known 
experimental values for total and partial decay widths of vector resonances 
to e+e" pair [Q. The calculation of background contribution at s > 1,44 



was performed as in the paper [0. The full integration region over s was 
divided on six intervals. In each interval the spectral density was founded 
by fitting of experimental points after subtraction of resonance contribution. 
Our numerical calculations show, that the main part of hadronic vacuum po- 
larization contribution to hyperfine splitting of muonic hydrogen caused by vr 
meson form factor (20). Theoretical error, not exceeding 10%, reflects spread 
in the choice of fitting parameters for cross section a{e^e~ — >■ hadrons) and 
for nucleon form factors. So, we actually see that whereas for hydrogen atom 
the relative contribution of hadronic vacuum polarization to ground state hy- 



perfine splitting is about 10 ^ [Q, for muonium - 5.6 -10 * [^, |T5[, then in the 
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case of muonic hydrogen it increases to 19.7 • 10^^. The essential modifica- 
tion of hadronic vacuum polarization relative contribution to HFS of muonic 
hydrogen is connected with increase of particle relative motion momentum, 
which achieves the values of order of hadronic masses. Obtained in this pa- 
per value of correction /S.E^^''^ = 861.05 MHz, must be accounted for future 
comparison with experimental result and for more exact determination of 6^ 
and 6^ [i. 

When constructing two-photon interaction quasipotential we take into ac- 
count in the virtual Compton scattering amplitudes only elastic proton inter- 
mediate states. We have investigated exactly this mechanism of the hadronic 
vacuum polarization contribution to HFS of muonic hydrogen, which gives 
the correction to 6^. It is possible additional "nonelastic" process consid- 
ering the excitement in the intermediate states of "ladder" diagrams of A 
isobar or other nucleon resonances (correction 6^). The contribution of such 
nonelastic mechanism demands special study. 



Interval y/s, Gev 


AE^^-' (fip), MHz 


P 


442.97 




45.74 





48.97 


J/^ family 


11.04 


T family 


0.11 


Background 




0.8<v^<1.2 R = A + Bs + Cs'^ 


131.90 


1.2 < < 3.0 R = 2.843s-°-isi 


104.01 


3.0 < < 7.4 R = 1.293s°-^^2 


70.78 


7.4 <v^< 12.0 R=3.91 


3.46 


12.0 < < 36.0 R=4.01 


1.94 


36.0 < < 47.0 R=4.13 


0.11 


^ > 47.0 QCD 


0.02 


Summary contribution 


861.05 



Nevertheless, it is possible to give the following approximate appreciation 
of hadronic vacuum polarization contribution to HFS of muonic hydrogen, 
taking into account proton polarizability. In the case of ordinary hydrogen 
the polarization correction in (2) is no more than 4 ppm [4]. It is equal 
approximately to 10% from elastic proton contribution 6^ . Then we can 
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obtain the approximate value of discussed contribution from proton polariz- 
ability replacing electron mass nie by muon mass (as it follows from (23)) 
and adding one degree of a for calculation of vacuum polarization. So, in the 
case of muonic hydrogen such " nonelastic" mechanism may give contribution 
to HFS of order of 4 ppm^a ^ 6.0 ppm. This value is about 30% from 
relative contribution of hadronic vacuum polarization, obtained in this paper 
and equal to 19.7ppm. 

We are grateful to F. Jegerlehner, A. Karimkhodzaev, S. Karshenboim 
and V. Saleev for useful discussions. This work was supported by the Russian 
Foundation for Basic Research (grant no. 96-02-17309). 
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